Single-cell genomics and regulatory networks for 388 human brains.

Single-cell genomics is a powerful tool for studying heterogeneous tissues such as the brain. Yet, little is understood about how genetic variants influence cell-level gene expression. Addressing this, we uniformly processed single-nuclei, multi-omics datasets into a resource comprising >2.8M nuclei from the prefrontal cortex across 388 individuals. For 28 cell types, we assessed population-level variation in expression and chromatin across gene families and drug targets. We identified >550K cell-type-specific regulatory elements and >1.4M single-cell expression-quantitative-trait loci, which we used to build cell-type regulatory and cell-to-cell communication networks. These networks manifest cellular changes in aging and neuropsychiatric disorders. We further constructed an integrative model accurately imputing single-cell expression and simulating perturbations; the model prioritized ~250 disease-risk genes and drug targets with associated cell types.

Disease-specific selective vulnerability and neuroimmune pathways in dementia revealed by single cell genomics.

The development of successful therapeutics for dementias requires an understanding of their shared and distinct molecular features in the human brain. We performed single-nuclear RNAseq and ATACseq in Alzheimer disease (AD), Frontotemporal degeneration (FTD), and Progressive Supranuclear Palsy (PSP), analyzing 40 participants, yielding over 1.4M cells from three brain regions ranging in vulnerability and pathological burden. We identify 35 shared disease-associated cell types and 14 that are disease-specific, replicating those previously identified in AD. Disease - specific cell states represent molecular features of disease-specific glial-immune mechanisms and neuronal vulnerability in each disorder, layer 4/5 intra-telencephalic neurons in AD, layer 2/3 intra-telencephalic neurons in FTD, and layer 5/6 near-projection neurons in PSP. We infer intrinsic disease-associated gene regulatory networks, which we empirically validate by chromatin footprinting. We find that causal genetic risk acts in specific neuronal and glial cells that differ across disorders, primarily non-neuronal cells in AD and specific neuronal subtypes in FTD and PSP. These data illustrate the heterogeneous spectrum of glial and neuronal composition and gene expression alterations in different dementias and identify new therapeutic targets by revealing shared and disease-specific cell states.

2-Deoxyglucose drives plasticity via an adaptive ER stress-ATF4 pathway and elicits stroke recovery and Alzheimer's resilience.

Intermittent fasting (IF) is a diet with salutary effects on cognitive aging, Alzheimer's disease (AD), and stroke. IF restricts a number of nutrient components, including glucose. 2-deoxyglucose (2-DG), a glucose analog, can be used to mimic glucose restriction. 2-DG induced transcription of the pro-plasticity factor, Bdnf, in the brain without ketosis. Accordingly, 2-DG enhanced memory in an AD model (5xFAD) and functional recovery in an ischemic stroke model. 2-DG increased Bdnf transcription via reduced N-linked glycosylation, consequent ER stress, and activity of ATF4 at an enhancer of the Bdnf gene, as well as other regulatory regions of plasticity/regeneration (e.g., Creb5, Cdc42bpa, Ppp3cc, and Atf3) genes. These findings demonstrate an unrecognized role for N-linked glycosylation as an adaptive sensor to reduced glucose availability. They further demonstrate that ER stress induced by 2-DG can, in the absence of ketosis, lead to the transcription of genes involved in plasticity and cognitive resilience as well as proteostasis.

Integrin-Driven Axon Regeneration in the Spinal Cord Activates a Distinctive CNS Regeneration Program.

The peripheral branch of sensory dorsal root ganglion (DRG) neurons regenerates readily after injury unlike their central branch in the spinal cord. However, extensive regeneration and reconnection of sensory axons in the spinal cord can be driven by the expression of α9 integrin and its activator kindlin-1 (α9k1), which enable axons to interact with tenascin-C. To elucidate the mechanisms and downstream pathways affected by activated integrin expression and central regeneration, we conducted transcriptomic analyses of adult male rat DRG sensory neurons transduced with α9k1, and controls, with and without axotomy of the central branch. Expression of α9k1 without the central axotomy led to upregulation of a known PNS regeneration program, including many genes associated with peripheral nerve regeneration. Coupling α9k1 treatment with dorsal root axotomy led to extensive central axonal regeneration. In addition to the program upregulated by α9k1 expression, regeneration in the spinal cord led to expression of a distinctive CNS regeneration program, including genes associated with ubiquitination, autophagy, endoplasmic reticulum (ER), trafficking, and signaling. Pharmacological inhibition of these processes blocked the regeneration of axons from DRGs and human iPSC-derived sensory neurons, validating their causal contributions to sensory regeneration. This CNS regeneration-associated program showed little correlation with either embryonic development or PNS regeneration programs. Potential transcriptional drivers of this CNS program coupled to regeneration include Mef2a, Runx3, E2f4, and Yy1. Signaling from integrins primes sensory neurons for regeneration, but their axon growth in the CNS is associated with an additional distinctive program that differs from that involved in PNS regeneration.SIGNIFICANCE STATEMENT Restoration of neurologic function after spinal cord injury has yet to be achieved in human patients. To accomplish this, severed nerve fibers must be made to regenerate. Reconstruction of nerve pathways has not been possible, but recently, a method for stimulating long-distance axon regeneration of sensory fibers in rodents has been developed. This research uses profiling of messenger RNAs in the regenerating sensory neurons to discover which mechanisms are activated. This study shows that the regenerating neurons initiate a novel CNS regeneration program which includes molecular transport, autophagy, ubiquitination, and modulation of the endoplasmic reticulum (ER). The study identifies mechanisms that neurons need to activate to regenerate their nerve fibers.

Molecular cascades and cell-type specific signatures in ASD revealed by single cell genomics.

Understanding how genetic variation exerts its effects on the human brain in health and disease has been greatly informed by functional genomic characterization. Studies over the last decade have demonstrated robust evidence of convergent transcriptional and epigenetic profiles in post-mortem cerebral cortex from individuals with Autism Spectrum Disorder (ASD). Here, we perform deep single nuclear (sn) RNAseq to elucidate changes in cell composition, cellular transcriptomes and putative candidate drivers associated with ASD, which we corroborate using snATAC-seq and spatial profiling. We find changes in cell state composition representing transitions from homeostatic to reactive profiles in microglia and astrocytes, a pattern extending to oligodendrocytes and blood brain barrier cells. We identify profound changes in differential expression involving thousands of genes across neuronal and glial subtypes, of which a substantial portion can be accounted for by specific transcription factor networks that are significantly enriched in common and rare genetic risk for ASD. These data, which are available as part of the PsychENCODE consortium, provide robust causal anchors and resultant molecular phenotypes for understanding ASD changes in human brain.

Transcription factor network analysis identifies REST/NRSF as an intrinsic regulator of CNS regeneration in mice.

The inability of neurons to regenerate long axons within the CNS is a major impediment to improving outcome after spinal cord injury, stroke, and other CNS insults. Recent advances have uncovered an intrinsic program that involves coordinate regulation by multiple transcription factors that can be manipulated to enhance growth in the peripheral nervous system. Here, we use a systems genomics approach to characterize regulatory relationships of regeneration-associated transcription factors, identifying RE1-Silencing Transcription Factor (REST; Neuron-Restrictive Silencer Factor, NRSF) as a predicted upstream suppressor of a pro-regenerative gene program associated with axon regeneration in the CNS. We validate our predictions using multiple paradigms, showing that mature mice bearing cell type-specific deletions of REST or expressing dominant-negative mutant REST show improved regeneration of the corticospinal tract and optic nerve after spinal cord injury and optic nerve crush, which is accompanied by upregulation of regeneration-associated genes in cortical motor neurons and retinal ganglion cells, respectively. These analyses identify a role for REST as an upstream suppressor of the intrinsic regenerative program in the CNS and demonstrate the utility of a systems biology approach involving integrative genomics and bio-informatics to prioritize hypotheses relevant to CNS repair.

Core transcription programs controlling injury-induced neurodegeneration of retinal ganglion cells.

Regulatory programs governing neuronal death and axon regeneration in neurodegenerative diseases remain poorly understood. In adult mice, optic nerve crush (ONC) injury by severing retinal ganglion cell (RGC) axons results in massive RGC death and regenerative failure. We performed an in vivo CRISPR-Cas9-based genome-wide screen of 1,893 transcription factors (TFs) to seek repressors of RGC survival and axon regeneration following ONC. In parallel, we profiled the epigenetic and transcriptional landscapes of injured RGCs by ATAC-seq and RNA-seq to identify injury-responsive TFs and their targets. These analyses converged on four TFs as critical survival regulators, of which ATF3/CHOP preferentially regulate pathways activated by cytokines and innate immunity and ATF4/C/EBPγ regulate pathways engaged by intrinsic neuronal stressors. Manipulation of these TFs protects RGCs in a glaucoma model. Our results reveal core transcription programs that transform an initial axonal insult into a degenerative process and suggest novel strategies for treating neurodegenerative diseases.

Toll-like receptor 2 (TLR2)-deficiency impairs male mouse recovery from a depression-like state.

Major depression is a prevalent, debilitating disease, yet therapeutic interventions for depression are frequently inadequate. Many clinical and pre-clinical studies have demonstrated that depression is associated with aberrant activation of the inflammatory system, raising the possibility that reducing inflammation may provide antidepressant effects. Using the learned helplessness mouse model, we tested if susceptibility or recovery were affected by deficiency in either of two receptors that initiate inflammatory signaling, Toll-like receptor-4 (TLR4) and TLR2, using knockout male mice. TLR4-/- mice displayed a strong resistance to learned helplessness, confirming that blocking inflammatory signaling through TLR4 provides robust protection against this depression-like behavior. Surprisingly, TLR2-/- mice displayed increased susceptibility to learned helplessness, indicating that TLR2-mediated signaling counteracts susceptibility. TLR2-mediated signaling also promotes recovery, as TLR2-/- mice demonstrated a severe impairment in recovery from learned helplessness. That TLR2 actually protects from learned helplessness was further verified by the finding that administration of the TLR2 agonist Pam3CSK4 reduced susceptibility to learned helplessness. Treatment with Pam3CSK4 also reversed chronic restraint stress-induced impaired sociability and impaired learning in the novel object recognition paradigm, demonstrating that TLR2 stimulation can protect from multiple impairments caused by stress. In summary, these results demonstrate that TLR2-mediated signaling provides a counter-signal to oppose deleterious effects of stress that may be related to depression, and indicate that TLR2 and TLR4 act oppositely to balance mood-relevant responses to stress.

Integrative methylome and transcriptome analysis to dissect key biological pathways for psoriasis in Chinese Han population.

BACKGROUND: Recent studies have revealed that DNA methylation (DNAm) could modulate gene expression in psoriasis (Ps). However, the relationship between whole-genome DNAm and gene expression in Ps has not been studied yet. OBJECTIVES: To better characterize the relationship between DNAm and gene expression, and to identify biological pathways triggered by changes in methylation involved in the pathogenesis of Ps. METHODS: Differentially methylated sites (DMSs) and differentially expressed genes (DEGs) were analysed by comparing 20 involved psoriatic (PP) skin, 20 uninvolved psoriatic (PN) skin and 20 normal (NN) skin biopsies. DEGs in negative correlation with the methylation were entered into further Gene Ontology (GO) and pathway analysis by clusterProfiler package in R program. RESULTS: A total of 290 genes with reverse correlation overlapped in PP vs PN and PP vs NN comparisons. GO categories of reversely-associated genes mainly enriched in T cell activation, type I interferon signaling pathway and defense response to other organism. Pathway analysis revealed superior NOD-like receptor signaling pathway and Measles enriched in the differentially up-regulated transcripts and regulation of lipolysis in adipocytes in the down-regulated transcripts. CONCLUSIONS: Our results provided a comprehensive correlation analysis of transcriptome and methylome in Ps. Increased innate immunity and decreased lipid biosynthesis play important roles in the development of psoriatic skin. This integrated analysis shed light on novel insights into the pathogenic mechanisms involved in Ps.

Novel mutations in Chinese Han patients with tuberous sclerosis complex: Case series and review of the published work.

Tuberous sclerosis complex (TSC) is an autosomal dominant genetic disease characterized by hamartomas in multiple organ systems. This study was performed in one familial and two sporadic cases with TSC. Two novel mutations (c.1884_1887delAAAG and c.5266A>G) and two previously reported mutations (c.4258_4261delTCAG and c.1960G>C) were identified by direct DNA sequencing. Of the four mutations, c.1884_1887delAAAG and c.1960G>C were found in a family and identified in the same allele by TA cloning sequencing. However, c.1960G>C was reported to be non-pathogenic. Furthermore, correlations between genotypes and phenotypes of Chinese Han patients since 2014 were performed by paired χ2 -tests in our published work review, which has not been reported. The results showed that patients with TSC2 mutations had a higher frequency of mental retardation and there were no significant differences of seizures and skin lesions with TSC1 mutations. Genetically, they had a higher frequency of familial inheritance.

Author Correction: Whole-exome SNP array identifies 15 new susceptibility loci for psoriasis.

This corrects the article DOI: 10.1038/ncomms7793.

TNFα disrupts blood brain barrier integrity to maintain prolonged depressive-like behavior in mice.

Recovery from major depressive disorder is difficult, particularly in patients who are refractory to antidepressant treatments. To examine factors that regulate recovery, we developed a prolonged learned helplessness depression model in mice. After the induction of learned helplessness, mice were separated into groups that recovered or did not recover within 4 weeks. Comparisons were made between groups in hippocampal proteins, inflammatory cytokines, and blood brain barrier (BBB) permeability. Compared with mice that recovered and control mice, non-recovered mice displaying prolonged learned helplessness had greater hippocampal activation of glycogen synthase kinase-3 (GSK3), higher levels of tumor necrosis factor-α (TNFα), interleukin-17A, and interleukin-23, increased permeability of the blood brain barrier (BBB), and lower levels of the BBB tight junction proteins occludin, ZO1, and claudin-5. Treatment with the GSK3 inhibitor TDZD-8 reduced inflammatory cytokine levels, increased tight junction protein levels, and reversed impaired recovery from learned helplessness, demonstrating that prolonged learned helplessness is reversible and is maintained by abnormally active GSK3. In non-recovered mice with prolonged learned helpless, stimulation of sphingosine 1-phosphate receptors by Fingolimod or administration of the TNFα inhibitor etanercept repaired the BBB and reversed impaired recovery from prolonged learned helplessness. Thus, disrupted BBB integrity mediated in part by TNFα contributes to blocking recovery from prolonged learned helplessness depression-like behavior. Overall, this report describes a new model of prolonged depression-like behavior and demonstrates that stress-induced GSK3 activation contributes to disruption of BBB integrity mediated by inflammation, particularly TNFα, which contributes to impaired recovery from prolonged learned helplessness.

Fine-mapping analysis of the MHC region for vitiligo based on a new Han-MHC reference panel.

Vitiligo is an immune-related disease with patchy depigmentation of skin and hair caused by selective destruction of melanocytes. In recent decades, many studies have shown the association between vitiligo and HLA genes; however, the results of Han Chinese are scarce. In this study, we performed a fine-mapping analysis of the MHC region in 2818 Han Chinese subjects through a widely used HLA imputation method with a newly built large-scale Han-MHC reference panel. Three new four-digit HLA alleles (HLA-DQB1 ∗ 02:02, HLA-DQA1 ∗ 02:01 and HLA-DPB1 ∗ 17:01) were identified to be associated with the risk of vitiligo, and four previously reported alleles were confirmed. Further conditional analysis revealed that two important variants, HLA-DQß1 amino acid position 135 (OR = 1.79, P = 1.87 × 10-11) and HLA-B amino acid positions 45-46 (OR = 1.44, P = 5.61 × 10-11), conferred most of the MHC associations. Three-dimension ribbon models showed that the former is located within the ß2 domain of the HLA-DQß1 molecule, and the latter lies in the α1 domain of the HLA-B molecule, while both are involved in specific antigen presenting process. Finally, we summarized all significant signals in the MHC region to clarify their complex relationships, and 8.60% of phenotypic variance could be explained based on all reported variants in Han Chinese so far. Our findings highlight the complex genetic architecture of the MHC region for vitiligo in Han Chinese population and expand our understanding of the roles of HLA coding variants in the etiology of vitiligo.

Exome-wide association study identifies four novel loci for systemic lupus erythematosus in Han Chinese population.

OBJECTIVES: Systemic lupus erythematosus (SLE) is a chronic autoimmune disease of considerable genetic predisposition. Genome-wide association studies have identified tens of common variants for SLE. However, the majority of them reside in non-coding sequences. The contributions of coding variants have not yet been systematically evaluated. METHODS: We performed a large-scale exome-wide study in 5004 SLE cases and 8179 healthy controls in a Han Chinese population using a custom exome array, and then genotyped 32 variants with suggestive evidence in an independent cohort of 13 246 samples. We further explored the regulatory effect of one novel non-coding single nucleotide polymorphism (SNP) in ex vivo experiments. RESULTS: We discovered four novel SLE gene regions (LCT, TPCN2, AHNAK2 and TNFRSF13B) encompassing three novel missense variants (XP_016859577.1:p.Asn1639Ser, XP_016859577.1:p.Val219Phe and XP_005267356.1:p.Thr4664Ala) and two non-coding variants (rs10750836 and rs4792801) with genome-wide significance (pmeta <5.00×10-8). These variants are enriched in several chromatin states of primary B cells. The novel intergenic variant rs10750836 exhibited an expression quantitative trait locus effect on the TPCN2 gene in immune cells. Clones containing this novel SNP exhibited gene promoter activity for TPCN2 (P=1.38×10-3) whose expression level was reduced significantly in patients with SLE (P<2.53×10-2) and was suggested to be further modulated by rs10750836 in CD19+ B cells (P=7.57×10-5) in ex vivo experiments. CONCLUSIONS: This study identified three novel coding variants and four new susceptibility gene regions for SLE. The results provide insights into the biological mechanism of SLE.

Genotype combination contributes to psoriasis: An exhaustive algorithm perspective.

Researchers have learned that nearly all conditions and diseases have a genetic component. With the benefit of technological advances, many single-nucleotide polymorphisms (SNPs) have been found to be associated with the risk of complex disorders by using genome wide association studies (GWASs). Disease-associated SNPs are sometimes shared by healthy controls and cannot clearly distinguish affected individuals from unaffected ones. The combined effects of multiple independent SNPs contribute to the disease process, but revealing the relationship between genotype and phenotype based on the combinations remains a great challenge. In this study, by considering the disease prevalence rate, we conducted an exhaustive process to identify whether a genotype combination pattern would have a decisive effect on complex disorders. Based on genotype data for 68 reported SNPs in 8,372 psoriasis patients and 8,510 healthy controls, we found that putative causal genotype combination patterns (CGCPs) were only present in psoriasis patients, not in healthy subjects. These results suggested that psoriasis might be contributed by combined genotypes, complementing the traditional modest susceptibility of a single variant in a single gene for a complex disease. This work is the first systematic study to analyze genotype combinations based on the reported susceptibility genes, considering each individual among the cases and controls from the Chinese population, and could potentially advance disease-gene mapping and precision medicine due to the causality relationship between the candidate CGCPs and complex diseases.

Novel Regulatory Mechanisms for the SoxC Transcriptional Network Required for Visual Pathway Development.

What pathways specify retinal ganglion cell (RGC) fate in the developing retina? Here we report on mechanisms by which a molecular pathway involving Sox4/Sox11 is required for RGC differentiation and for optic nerve formation in mice in vivo, and is sufficient to differentiate human induced pluripotent stem cells into electrophysiologically active RGCs. These data place Sox4 downstream of RE1 silencing transcription factor in regulating RGC fate, and further describe a newly identified, Sox4-regulated site for post-translational modification with small ubiquitin-related modifier (SUMOylation) in Sox11, which suppresses Sox11's nuclear localization and its ability to promote RGC differentiation, providing a mechanism for the SoxC familial compensation observed here and elsewhere in the nervous system. These data define novel regulatory mechanisms for this SoxC molecular network, and suggest pro-RGC molecular approaches for cell replacement-based therapies for glaucoma and other optic neuropathies.SIGNIFICANCE STATEMENT Glaucoma is the most common cause of blindness worldwide and, along with other optic neuropathies, is characterized by loss of retinal ganglion cells (RGCs). Unfortunately, vision and RGC loss are irreversible, and lead to bilateral blindness in ∼14% of all diagnosed patients. Differentiated and transplanted RGC-like cells derived from stem cells have the potential to replace neurons that have already been lost and thereby to restore visual function. These data uncover new mechanisms of retinal progenitor cell (RPC)-to-RGC and human stem cell-to-RGC fate specification, and take a significant step toward understanding neuronal and retinal development and ultimately cell-transplant therapy.

Intranasal siRNA administration reveals IGF2 deficiency contributes to impaired cognition in Fragile X syndrome mice.

Molecular mechanisms underlying learning and memory remain imprecisely understood, and restorative interventions are lacking. We report that intranasal administration of siRNAs can be used to identify targets important in cognitive processes and to improve genetically impaired learning and memory. In mice modeling the intellectual deficiency of Fragile X syndrome, intranasally administered siRNA targeting glycogen synthase kinase-3ß (GSK3ß), histone deacetylase-1 (HDAC1), HDAC2, or HDAC3 diminished cognitive impairments. In WT mice, intranasally administered brain-derived neurotrophic factor (BDNF) siRNA or HDAC4 siRNA impaired learning and memory, which was partially due to reduced insulin-like growth factor-2 (IGF2) levels because the BDNF siRNA- or HDAC4 siRNA-induced cognitive impairments were ameliorated by intranasal IGF2 administration. In Fmr1-/- mice, hippocampal IGF2 was deficient, and learning and memory impairments were ameliorated by IGF2 intranasal administration. Therefore intranasal siRNA administration is an effective means to identify mechanisms regulating cognition and to modulate therapeutic targets.